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ABSTRACT 

The  Delaware  Method  of  rating  baffled  shell  and 
tube  heat  exchangers  is  used  as  the  basis  for  a  simplified 
mathematical  model  for  predicting  shell  side  heat  transfer 
coefficients  and  pressure  drops.   Multiple  non-linear 
regression  analysis  is  used  to  express  the  heat  transfer 
coefficient  and  pressure  drop  as  a  product  of  terms 
involving  the  five  basic  geometric  parameters  of  a  shell- 
and-tube  heat  exchanger,  over  a  suitable  range  of  these 
variables. 

A  brief  sensitivity  analysis  of  heat  transfer  per 
unit  pressure  drop  is  presented,  along  with  a  proposed 
routine  for  determining  exchanger  size  when  the  shell- 
side  coefficient  is  the  controlling  factor.   A  comparison 
is  made  between  this  method  and  the  Delaware  Method. 
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NOMENCLATURE 

A.  heat  transfer  area  for  tube  side,  (ft  ) 

A  heat  transfer  area  for  shell  side,  (ft  ) 
o 

c  specific  heat  of  shell  side  fluid,  (Btu/lbm  -  °F) 

D-  shell  inside  diameter,  (in) 

D  . „  diameter  of  outer  limit  of  tube  bundle,  (in) 

d  tube  outside  diameter,  (in) 

d.  tube  inside  diameter,  (in) 

F  fraction  of  tubes  that  are  in  crossflow 
c 

F  .  fraction  of  total  crossflow  area,  available  for 

s  "  bypass  flow  around  tube  bundle 


f.      gravitational  friction  factor  for  flow  across 
an  ideal  tube  bank 

o 

g       gravitational  conversion  constant,  4.17  x  10 
c      lbm-ft/lbf-hr2 

h.      tube  side  heat  transfer  coefficient  (Btu/hr-ft  -  F) 

h-j   ,   ideal  shell  side  heat  transfer  coefficient, 
ldeal   (Btu/hr-ft2-0F) 


h       exchanger  shell  side  heat  transfer  coefficient, 
0      (Btu/hr-ft2-°F) 
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J,       shell  side  heat  transfer  coefficient  correction 
factor  for  bundle  bypass  effects 


J       shell  side  heat  transfer  coefficient  correction 
factor  for  baffle  configuration  effects 


J.      shell  side  heat  transfer  coefficient  correction 
factor  for  baffle  leakage  effects 


J       shell  side  heat  transfer  coefficient  correction 
factor  for  adverse  temperature  gradients 


J       shell  side  heat  transfer  coefficient  correction 
factor  for  the  effect  of  unequal  baffle  spacing 


j.      Colburn  j -factor  for  ideal  tube  bank 

K,       thermal  conductivity  of  tube  side  fluid, 
D      (Btu/hr-ft-°F) 


K       thermal  conductivity  of  shell  side  fluid, 
s      (Btu/hr-ft-°F) 

K       thermal  conductivity  of  tube  wall,  (Btu/hr-ft-  F) 


w 


.5T.     logarithmic  mean  temperature  difference 

L       effective  tube  length  between  tube  sheets,  (ft) 

I  baffle  cut  distance  from  baffle  to  shell  inside 


c 


diameter,  (in) 


I  baffle  spacing,  (in) 


I      ■ , i  baffle  spacing  at  inlet  and  outlet  of  exchanger, 

s,i   s,   respectively 


N,        number  of  baffles  in  exchanger 

N         number  of  tube  rows  crossed  during  flow  through 
one  crossflow  section 


N        number  of  effective  crossflow  rows  in  each 
cw       window  section 


N         total  number  of  tubes  in  exchanger 

APv  •     pressure  drop  for  flow  across  one  ideal  crossflow 
0,1     section,  (lbf/in2) 


AP   •     pressure  drop  for  flow  through  one  ideal  window 
w>1     section,  (lbf/in2) 


p       •  tube  pitch,  (in) 


p         tube  pitch  normal  to  flow,  (in) 

p         tube  pitch  parallel  to  flow,  (in) 

Q  '  total  heat  transferred  in  heat  exchanger, 

(Btu/hr) 


R,        pressure  drop  correction  factor  for  effect  of 
bundle  bypass 


R         pressure  drop  correction  factor  for  effect  of 
baffle  leakage 


•.&' 


R         pressure  drop  correction  factor  for  effect  of 


s 


unequal  inlet  and  outlet  baffle  spacing 
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R         shell  side  Reynolds  number 
es 

S         crossflow  area  at  centerline  for  one  crossflow 


m 


o 


section,  (in^) 


S  .       shell-to-baffle  leakage  area  for  one  baffle,  (in  ) 

S  ,       tube-to-baffle  leakage  area  for  one  baffle,  (in  ) 

2 
S         area  for  flow  through  window,  (in  ) 
w  e 

2 
S         total  baffle  window  area,  (in  j 
wg 

2 
S  .       window  area  occupied  by  tubes,  (in  ) 

U         overall  heat  transfer  coefficient  base  on  shell 


side  heat  transfer  area,  (Btu/hr-ft^-OF) 


W  mass  flow  rate  of  shell  side  fluid,  (lbm/hr) 

W  mass  flow  rate  of  tube  side  fluid,  (lbm/hr) 

<5  ,  diametral  clearance  between  shell  and  baffle, 

Sb  (in) 


6.,       diametral  clearance  between  tube  and  baffle, 
tb       (in) 


6   ,      approximate  diametral  clearance  between  shell 
and  outer  tube  limit,  (in) 


y         shell  side  fluid  viscosity  at  bulk  temperature, 


s 


(lbm/hr-ft) 


-10 


y         shell  side  fluid  viscosity  at  surface  temperature, 
s'w      (lbm/hr-ft) 

3 
p         shell  side  fluid  density,  (lbm/ft  ) 

9         baffle  cut  angle,  (rad) 


P  shell  side  Prandtl  number 
rs 

y,  tube  side  fluid  viscosity  at  bulk  temperature, 

D  (lbm/hr-ft) 
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I.   INTRODUCTION 

A.       BACKGROUND 

The  aim  of  any  heat  exchanger  rating  method  is  to 
predict,  as  accurately  as  possible,  the  performance  of  a 
heat  exchanger  with  respect  to  total  heat  transferred  and 
pressure  drops  encountered  on  both  the  shell  and  tube  sides 

The  Delaware  Method  for  Rating  Baffled  Shell  and 
Tube  Heat  Exchangers  (1)  was  developed  from  its  inception 
in  the  late  1940 's  through  the  early  1960 's  as  a  result  of 
numerous  experiments  conducted  at  the  University  of 
Delaware.   This  method  is  applicable  to  many  variations 
on  the  baffled  exchanger:   multiple  passes,  finned  tubes, 
various  tube  orientations,  a  wide  range  of  Reynolds 
numbers,  etc.   As  Bell  has  stated   (1),  the  method  does 
have  its  uncertainties,  as  all  methods  do.   However,  this 
method  is  widely  used  and  easy  to  understand. 

Tinker  (2)  and  Devore  (3)  propose  similar  methods 
for  rating  baffled  exchangers,  but,  in  the  opinion  of 
this  author,  are  not  as  straightforward  as  the  Delaware 
Method.   Bell  himself  has  proposed  routines  for  sizing 
baffled  exchangers  based  on  the  Delaware  Method 
(references  1  and  11) ,  and  this  method  has  been  used 
extensively  in  computer  programs  for  "optimal"  exchanger 
determination.   References (9)  and  (10)  are  two  such 
examples . 
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The  basic  means  of  analysis  in  the  Delaware  Method 
takes  into  consideration  five  fluid  streams  and  three 
heat  transfer  zones  in  the  baffled  exchanger,  as  illustrated 
in  Figures  1  and  2.   As  can  be  seen  by  these  simplified 
illustrations,  the  flow  pattern  for  shell  side  fluid  is 
complex  at  best.   Reference  (1)  contains  the  complete  set 
of  equations  for  relating  the  heat  transfer  coefficient 
and  pressure  drops  to  the  geometry  of  the  exchanger.   This 
reference  also  utilizes  extensively  graphs  and  tables 
for  determination  of  some  parameters.   Mueller  (7)  has 
developed  equations  for  many  of  these  graphical  relations. 
Appendix  A  is  a  compilation  of  the  basic . equations  of  the 
Delaware  Method  and  those  derived  by  Mueller.   This 
compilation  does  have  its  limitations  as  noted  in  Section  I.C 


FIGURE  1   Schematic  Diagram  of  Streams  on 
Shell  Side  of  Baffled  Heat 
Exchangers . 
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FIGURE  2(a)    Region  o£  Crossflow  Between  Baffle 
Tips  in.  the  Central  Baffle  Spacing 
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FIGURE  2(b)    Flow  Region  Considered  for  Window  Flow 


FIGURE  2(c)    Flow  Region  for  the  End  Baffle  Spacings 
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B.       PURPOSE  AND  METHODOLOGY 

The  purpose  of  this  study  is  to  provide  an  insight 
into  the  relations  between  shell-side  heat  transfer 
coefficient  and  pressure   drop  and  the  five  basic 
exchanger  parameters:   Tube  diameter  --  d  •  pitch-to- 
diameter  ratio  --  P/d  :  baffle  cut  ratio  --  £„/D.  ;  baffle 

o  c   i ' 

spacing  ratio  --  £  /D-;  and  shell  diameter  --  D.. 

It  is  anticipated  that  a  simplified  model  for 
determining  the  heat  transfer  coefficient  will  be  of 
considerable  use  to  the  designer  as  a  "first  cut"  in 
sizing  an  exchanger  for  a  given  job.   Admittedly,  this 
cannot  replace  the  sophisticated  computer  programs 
available  for  the  determination  of  "optimal"  exchanger 
configuration  for  a  given  application.   However,  the 
method  developed  here  will  be  a  step  in  the  right 
direction  for  initial  selection  of  parameters  and  will 
provide  some  "feel"  for  the  effect  of  geometry  variations 
on  overall  heat  exchanger  performance. 

This  analysis  was  conducted  utilizing  the  basic 
equations  of  the  Delaware  Method  as  modified  by  the 
assumptions  and  approximations  enumerated  in  Section  I.C. 
Utilizing  these  equations,  those  factors  affecting  heat 
transfer  and  pressure  drop  which  are  a  function  of 
exchanger  geometry  were  calculated  over  a  suitable  range 
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of  variables.   The  resulting  values  were  correlated 
utilizing  four  FORTRAN  subroutines  (b) :   CORRE,  ORDER, 
MINV  and  MULTR.   The  resulting  non-linear  correlation 
yields  these  factors  expressed  as  a  product  of  the 
five  basic  exchanger  parameters,  each  with  its  own 
regression  coefficient  as  an  exponent,  and  a  single 
correlation  coefficient.   These  regressions  provided 
acceptable  correlation  results,  the  maximum  overall 
error  being  on  the  order  of  five  percent. 

C.       LIMITATIONS/ASSUMPTIONS/APPROXIMATIONS 

1.    Limitations 

In  this  analysis,  no  attempt  has  been  made  to 
correlate  data  for  the  shell  side  heat  transfer  co- 
efficient where  sealing  strips  are  used  to  limit  the 
bundle  bypass  flow.   Nor  has  an  attempt  been  made  with 
regard  to  enhanced  surfaces  for  heat  transfer,  i.e. 
finned  tubes.   Finally,  only  the  single  pass  condition 
for  the  tube  side  fluid  is  taken  into  consideration. 
Admittedly,  these  restrictions  limit  the  scope  of 
this  study  with  regard  to  overall  applicability  to 
heat  exchanger  design.   However,  it  is  hoped  that  this 
analysis  will  provide  insight  into  the  basic  relations 
between  shell  side  heat  transfer  coefficient  and  the 
fundamental  exchanger  goemetric  parameters. 


-17- 


2.    Assumptions 

It  is  realized  that  manufacturing  considerations 
effect  significantly  the  design,  construction  and  final 
form  of  any  heat  exchanger.   Reference  (5)  sets  standards 
for  the  tube-to-baffle  and  baffle-to-shell  diametral 
clearances  for  various  applications.   These  values  are 
predominantly  of  the  "step  function"  type,  in  that  they 
can  not  be  expressed  as  a  continuous  dependent  variable 
based  on  any  of  the  independent  exchanger  parameters. 
For  the  purposes  of  this  study,  both  clearances  were 
assumed  constant  throughout  the  range  of  parameters 
considered.   The  values  used  for  these  clearances  are: 

Shell-to-baffle  clearance:   6  ,  =  0.175  inches 

sb 

Tube-to-baffle  clearance:    5  .  =  0.03125  inches 

These  assumptions  again  allow  concentration  on  the  effect 
of  major  exchanger  parameters. 

The  second  major  assumption  is  the  neglect  of 
unequal  baffle  spacing  for  inlet  and  outlet  passages. 
This  assumption  is  made  in  both  heat  transfer  coefficient 
and  pressure  drop  calculations.   Unequal  baffle  spacing, 
in  an  exchanger  where  the  ratios 


lii  =  -111   =  1.5, 


s      s 
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contributes  to  an  approximate  2-3%  reduction  in  the  over 
all  heat  transfer  coefficient  for  a  sufficiently  large 
number  of  baffles.   Where  pressure  drop  is  concerned, 
the  assumption  must  be  made  with  the  realization  that 
the  end  zones  be  considered  as  crossflow  areas  identical 
to  the  zones  between  the  baffles.   Equations  (19)  and 
(25)  of  Appendix  A  reflect  these  changes. 

3.    Approximations 

In  order  to  implement  any  numerical  analysis  of 
the  type  used  in  this  study,  it  is  necessary  to  derive 
relations  from  data  provided  by  the  manufacturing 
community  and  from  curves  of  friction  factor  and 
Colburn  j -factor. 

For  the  former,  references  (1),  (5)  and  (8)  were 
used  for  two  such  derivations.   The  number  of  tubes  in  a 
heat  exchanger  of  particular  shell  diameter  was  corre- 
lated as  a  function  of  shell  inside  diameter,  pitch- 
to-tube  diameter,  and  tube  diameter,  for  the  single- 
pass  condition.   These  results  are: 

For  square/rotated  square  tube  layouts: 


N  =  0.3559 


D. 

l 


2.17526 


(d   1 
o 


CD 


-19- 


For  triangular  tube  layouts 


N  =  0.40759 


D. 

1 


2.1797 


^0^ 


(2) 


These  approximations  give  results  within  5%    over  the  range 
of  shell  diameters  from  10  to  60  inches.   The  clearance 
for  bundle  bypass  flow  (D-  -  D ot»)  was  correlated  as  a 
function  of  shell  diameter  also,  resulting  in  the 
approximation: 


6  «,  =  (D.  -  D  .„)  =  0.5817  D. 
stA    *■  i    otlJ  l 


0.3066 


(3) 


Figures  3,  4,  and  5  are  graphs  of  friction  factor 
and  Colburn  j -factor  versus  shell-side  Reynolds  number 
for  the  various  tube  orientations  at  varying  pitch-to- 
tube  diameter  ratios.   Correlation  results  are  presented 
in  Tables  1  and  2  for  friction  factor  and  Colburn  j -factor, 
respectively,  for  the  separate  ranges  of  Reynolds  numbers 
shown. 

Table  3  presents  relations  between  tube  pitch, 
tube  pitch  normal  to  fluid  flow  and  tube  pitch  parallel 
to  fluid  flow  for  the  various  orientations.   These 
relations  are  used  in  equations  (1),  (3),  and  (4)  of 
Appendix  A  in  calculating  the  various  factors  associated 
with  heat  transfer  and  pressure  drop  calculations. 
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P  =  tube  pitch 

P  =  tube  pitch  normal  to  flow 

P  =  tube  pitch  parallel  to  flow 
P 


Tube  p  p 

Orientation  n  p 


square  1.0  P  1.0  P 


rotated  ^       p  ^ 

square  ' 


triangular  P/2  /T/2  P 


TABLE  3    Relations  Between  Tube  Pitch 
P,  and  P-, Pp  for  Various  Tube 
Orientations 
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D.      NEGLECTED  FACTORS  Jr  and  Js 

For  preliminary  estimation  and  sizing  as  proposed 
by  this  simplification,  the  factors  J   (correction  due  to 
adverse  temperature  gradients)  and  J   (correction  due  to 
unequal  inlet  and  outlet  baffle  spacing)  have  been 
neglected. 

To  account  for  the  effect  of  J  ,  the  final  heat 
exchanger  length  must  be  known,  which  allows  determination 
of  the  end  baffle  spacing.   J  may  then  be  calculated 


from: 


_  v  b-1^ s  ,  1 s  ,  o 


where 


j   =  _±L^ ^ _£_^ (4) 

^Nb-1)  +  *s,i  +  S,0 


*s,i  ~  ls,i/ls>        ls>o        ls,o/Zs 


The  initial  value  of  h  must  then  be  multiplied  by  J   to 
give  a  more  accurate  estimate. 

The  factor  J   is  calculated  only  in  laminar  flow 
(Re   <  100)  from: 

1.51 

J   =  0.18  (5) 

(Nc) 
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Should  this  condition  (Re   <  100)  exist,  again  adjust  h 
by  multiplying  by  J  .   Figures  6  and  7  show  graphically 
the  relations  for  J  vs  I      ■    and  J  vs  N  .   These  graphs 
may  be  used  in  lieu  of  actual  calculations. 
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FIGURE  6    Correction  Factor  for  Unequal  Baffle  Spacing, 

J   vs  I*    .  (  =  £*  J 
S      S  ,1  v   s  ,oy 
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1.51 


(Nc) 


0.18 
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LO- 


0.95-- 


O.9.. 


0.85- 


10 


15 


20 


25 


N 


FIGURE  7    Correction  Factor  for  Adverse  Temperature 

Gradients,  J   vs  N 
'   r     c 


so- 


il.   CORRECTION  FACTORS  INFLUENCING  THE 
HEAT  TRANSFER  COEFFICIENT 

A.       BUNDLE  BYPASS  FACTOR  --  Jb 

The  bundle  bypass  factor,  J,  ,  accounts  for  the 
effect  of  fluid  flow  between  the  shell  and  outer  tubes 
of  the  tube  bundle.   This  portion  of  the  flow  bypasses 
the  main  crossflow  section  of  the  tube  bank  and  has  a 
significant  effect  on  heat  transfer.   Equations  A-4, 
A-5,  and  A-18  are  used  to  calculate  Jr.  where  the  shell- 
to-tube-bundle  clearance  (D-  -  D  t,)  has  been  approxi- 
mated by  equation  (3).   Equation  A-18  was  developed  by 
Mueller  (7)  where  for  laminar/transition  flow  C,,  =  1.35 
and  for  turbulent  flow  C,,  =  1.25.   For  purposes  of  this 


considered:   1  <  RP   <  10   for  laminar/transition  flow 

cs 


analysis  only  two  regimes  of  Reynolds  numbers  are 

msidered:   1  <  RP   <  10   for  laminar/ 

cs 

and  103  <  Re   <  105  for  turbulent  flow. 
cs 

The  non-linear  regression  results  for  J,  are 

presented  in  Table  4  and  a  representative  curve  is  shown 

in  Figure  8.   As  can  be  seen,  J,  can  be  represented  as 

a  function  of  shell  diameter  D-,  tube  diameter  d  ,  and 

1  o 

pitch-to-diameter  ratio  P/d  . 

From  Table  4,  it  can  be  readily  seen  that  tube 

diameter  has  little  effect  on  J,  .   Over  the  range  of 

b  s 

d  considered  (0.5<d  <1.5),  the   change  in  J,  due  to 
o  o  b 
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increasing  (or  decreasing)  d   from  1.0  to  1.5  (or  to  0.5) 
is  less  than  1%  . 

Figure  8  is  applicable  only  for  turbulent  flow 
over  square  tube  orientations.   However,  we  may  estimate 
J,  for  other  flow/orientation  conditions  easily  from 
Table  4  and  can  estimate  changes  in  J,  in  the  following 
fashion.   If  J,  is  desired  for  rotated-square  layouts 
in  turbulent  flow,  we  may  form  a  ratio: 


a      a 

Jb    [C   (P/d  )  2  (D.)  3]r 
r     o  _    oJ         i    iJ 

a      a 
Jbs    CC0  (P/d0)  2  (Dt)  3]S 


(P/d0) 


»  A  3-  _ 

2        s 

s  CD±D  r 


(6) 


This  simple  relation  yields  the  result  that  rotated-square 
layouts,  for  P/d  =  1.5,  D.  =  30,  result  in  a  bundle 
bypass  factor  approximately  t%   higher  than  square  tube 
orientations. 

Such  insight  into  the  intricacies  of  the  Delaware 
Method  can  allow  the  designer  to  make  reasonable  "fine- 
tuning"  corrections  to  a  preliminary  design  without  the 
need  for  costly  (and  sometimes  cumbersome)  computer  models 
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a-        a2      a_ 

J   =  c   (d  )  l    (P/d  )  z  (D.) 
b     o    cr    v    o       1 


R 
es 

Layout 

c 

0 

al 

a2 

a, 

0 

Square 

.335 

-.0137 

.5249 

.1806 

1  -  103 

Rotated- 
Square 

.416 

-.0119 

.4358 

.1493 

Triang. 

.496 

-.0099 

.3503 

.1194 

Square 

.369 

-.0125 

.4789 

.1648 

3     5 

io°  -  lcr 

Rotated- 
Square 

.445 

-.0108 

.3976 

.1362 

Triang. 

.527 

-.009 

.3197 

.1089 

TABLE  4   Correlation  Results  for  Baffle 
Configuration  Factor,  J, 
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FIGURE  8    Bundle  Bypass  Factor  Jfe   vs   P/dQ 
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B.       BAFFLE  CONFIGURATION  FACTOR  --  Jc 

Equations  A-2  and  A-16  represent  the  baffle  con- 
figuration factor,  J  ,  which  amounts  for  changes  in  heat 
transfer  due  to  those  tubes  physically  located  in  the 
baffle  window  and  not  subjected  to  full  corssflow  from 
the  main  stream  of  fluid.   Equation  A-16  is  a  linear 
approximation  to  the  graph  of  J   developed  by  Bell.   Using 
equation  (3]  for  D  t-,  J   can  be  reasonably  approximated 
as  a  product  of  shell  diameter  and  baffle  cut  (or  baffle 
cut  ratio) .   This  result  is  independent  of  flow  condition 
and  tube  orientation.   Regression  analysis  results  in 
the  following  equation  for  J  : 

-0.5467      -0.0184 
Jc  =  0.492  (£C/Di)  (Di)  (7) 


Figure  9  is  a  family  of  curves  plotting  &r/D-  vs  J   at 

various  values  of  D. . 

i 

As  is  easily  seen  from  this  figure,  the  dependence 
of  J  on  shell  diameter  is  minimal  when  compared  with 
that  of  the  baffle  cut  ratio.   For  the  range  of  baffle 
cut  ratios  considered  in  this  analysis  (0.2  <  £  /D.  <  0.4), 
the  heat  transfer  coefficient  may  vary  from  a  15%  increase 
to  a  25%  decrease.   Bell  (1)  states  that  for  "well-designed" 
exchangers,  J   should  be  about  1.0,  which  corresponds  to 
a  value  of  &-/D-  of  about  0.25. 
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FIGURE  9    Baffle  Configuration  Factor  Jc  vs  Baffle 
Cut  Ratio  ^C/Di 
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It  is  interesting  to  note  that  in  the  optimization 
routines  of  references  (9)  and  (10),  the  baffle  cut  ratios 
for  "optimal"  design  were  0.25  and  0.2,  respectively. 
Admittedly,  the  objective  functions  for  these  optimizations 
were  slightly  different  (heat  transfer  area  and  total 
tube  length,  resp.),  but  the  similarity  of  results  is 
significant. 

C.       BAFFLE  LEAKAGE  FACTOR  --  J^ 

The  baffle  leakage  factor,  J-,  is  by  far  the  most 
complex  of  all  the  factors  affecting  the  overall  heat 
transfer  coefficient.   It  accounts  for  the  effect  of 
leakage  through  the  baffle  due  to  clearances  between 
baffle  and  shell  and  between  tube  and  baffle.   In  this 
regard,  it  is  a  function  of  all  five  basic  exchanger 
parameters.   Equations  A-6,  7,  8,  15,  16,  17  must  all 
be  used  in  the  calculation  of  J„,  where  some  additional 
dependence  is  noted  as  a  function  of  tube  orientation. 

The  non-linear  regression  results  for  J«  are 
presented  in  Table  5.   Numerous  graphs  and  families  of 
curves  can  be  generated  using  these  results.   However, 
for  the  sake  of  simplicity,  Figures  10a,  b,  and  c  are 
represenative  of  the  effect  on  J   of  several  parameter 
variations.   Bell  states  that  the  value  of  J.  should 
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3.  3.  3.  3. 

h   =  W  1  <P/V  2  (W  3  <VDi>  *  <V 


a5 


Layout      c       a,       a2       a.       a.       a^ 


Square     0.572    0.1296    0.7745    0.0975    0.2708    0.0718 
0.613    0.1177    0.7136    0.0856    0.2339    0.0643 
Triangular   0.656    0.1127    0.6679    0.0767    0.2051    0.0516 


Rotated- 
Square 


TABLE  5   Correlation  Results  for  Baffle  Leakage 
Factor,  J. 
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FIGURE  10a   Baffle  Leakage  Factor  J   vs  P/d 
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FIGURE  10b    Baffle  Leakage  Factor  J £   vs   P/dQ 
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FIGURE  10c    Baffle  Leakage  Factor  J£   vs   P/dQ 
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be  0.7  --  0.8  for  well-designed  exchangers.   From 
Figure  10,  one  can  easily  read  an  appropriate  value 
of  the  necessary  parameter  that  will  be  a  reasonable 
starting  point  for  heat  exchanger  design.   For  other 
shell  diameters  and  tube  diameters,  the  relations  of 
Table  5  may  be  suitably  manipulated  to  arrive  at 
reasonable  values  of  the  necessary  parameters. 

D.       MULTIPLE  CORRELATION  RESULTS 

To  facilitate  the  overall  prediction  of  the 

heat  transfer  coefficient  as  a  function  of  those 

variables  considered,  the  actual  shell  side  coefficient, 

h  ,  was  calculated  as  a  product  of  h. ,   ,  and  the 
o'  r  ideal 

aforementioned  factors:  J,  ,  J  ,  J„.   Colburn  i -factor 

b   c   I  J 

correlations  of  Table  1  were  used  for  the  various  flow 

regimes  as  were  correlations  for  H.    and  <5  .  „  .   The 
°  t      stx 

fluid  properties  and  mass  flow  rate  terms  were  factored 
out  of  the  equation  and  only  those  terms  involving 
geometric  parameters  were  used  for  the  overall  cor- 
relation.  These  terms  were  then  correlated  using  the 
aforementioned  Fortran  subroutines.   The  correlation 
results  are  presented  in  Table  6  for  the  various  flow 
regimes  and  tube  orientations.   It  should  be  noted  that 
all  input  dimensions  for  the  exchanger  parameters  are 
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in  inches  and  that  the  correction  factors  inherent  in 
these  equations  have  been  accounted  for  in  the  cor- 
relation constant.   The  resulting  units  for  h  are 
those  typically  expected  of  the  English  system,  Btu/HR- 
FT  -  F.   The  units  for  fluid  properties  are  as  listed 
in  the  nomenclature. 

Table  6  is  a  "busy"  listing  and  many  graphical 
representations  can  be  drawn  plotting  the  heat  transfer 
coefficient  h  against  the  various  parameters.   That 
will  not  be  done  here  but  some  aspects  of  Table  6  are 
noteworthy. 

Most  exchangers  are  designed  for  operation  in 

turbulent  flow  and  for  this  regime  it  can  be  seen  that 

h  is  nearly  inversely  proportional  to  shell  diameter, 

the  exponent  a,,  being  close  to  -1.0.   This  is  a  striking 

result,  in  that  for  desired  increases  in  h  ,  changes  in 
7  o 

shell  diameter  maybe  the  most  rapid  means  of  attaining 

the  desired  value. 

The  dependence  of  h  on  baffle  cut  ratio, 
r  o  ' 

I   /D. ,  is  independent  of  flow  but  dependent  on  tube 

orientation,  and  h  varies , nearly  as  the  square-root 

of  I   /D..   Pitch-to-diameter  ratio,  P/d  ,  has  little 
c   l  '    o* 

effect  on  square  tube  orientations,  but  can  cause 

significant  changes  in  triangular  and  rotated  square 

patterns.   The  dependence  of  h  on  tube  diameter  and 

r  o 
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baffle  spacing  ratio  is  straightforward  and  needs  no 
further  explanation. 

From  Table  6  it  is  quite  easy  to  predict  the 
heat  transfer  coefficient  given  the  exchanger  geometric 
parameters  and  fluid  properties.   Only  a  quick  Reynolds 
number  calculation  is  necessary  beforehand  to  determine 
the  applicable  flow  regime.   Changes  in  h   are  pre- 
dicted simply  by  calculating  the  ratio  of  the  parameter 
varied  raised  to  the  appropriate  exponent.   This  is 
quickly  accomplished  if  only  one  parameter  is  varied 
at  a  given  time.   For  large  changes  one  must  be  cautious 
that  the  flow  regime  remains  as  desired,  i.e.,  turbulent 
or  laminar/transition. 
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III.   CORRECTIONS  AFFECTING  PRESSURE  DROP 

A.       BAFFLE  LEAKAGE  FACTOR  --  R_£ 

The  pressure  drop  correction  factor  for  baffle 
leakage,  R^,  accounts  for  pressure  drop  changes  in  much 
the  same  fashion  that  J-  accounts  for  changes  in  heat 
transfer.   It  relates  the  leakage  through  the  baffle  due 
to  tube-to-baffle  and  baffle-to-shell  clearances  to  the 
geometric  parameters  of  the  heat  exchanger.   Equations 
A-2,  4,  6,  7,  8,  14,  15,  20  and  21  are  all  used  to 
calculate  R  .   Dependence  on  tube  orientation  is  accounted 
for  in  equations  A-4  and  A-6  by  the  use  of  the  correla- 
tions for  N   in  Equations  (1)  and  (2)  and  the  tube  pitch 
relations  of  Table  3. 

The  non-linear  correlation  results  are  presented 
in  Table  7  and  a  representative  graph  as  Figure  11. 
Figures  10  and  11  show  a  striking  resemblance  of  form, 
which  is  to  be  expected.   One  can  logically  assume  that 
an  increase  in  heat  transfer  coefficient  due  to  less 
leakage  through  the  baffle  will  also  result  in  increased 
pressure  drops. 
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a  3.  a  a 


a5 


Layout       aQ      &1      a£       a3      a4      a5 


Square  .2985  .2217  1.315  .1754  .4674  .1584 

Rotated     32g  >17gg  1>147  >1475  .3954  %157 
Square 

Triangular  .3702  .1682  1.0696  .1318  .3496  .1365 


TABLE  7   Correlation  Results  for  Baffle  Leakage  Factor  R. 
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B.       BUNDLE  BYPASS  FACTOR  --  Rb 

The  bundle  bypass  factor,  R,  ,  accounts  for  pressure 
drop  corrections  in  the  same  manner  that  J,  accounts  for 
heat  transfer  changes.   From  Equation  A-22,  it  can  be 
seen  to  relate  that  correction  to  F  ,  ,  the  fractional 
crossflow  area  available  for  bundle  bypass. 

Table  8  presents  the  correlation  results  for  R,  . 
Similar  to  J,  is  the  weak  dependence  on  tube  diameter 
and  the  strong  dependence  on  pitch-to-diameter  ratio. 
This  comes  as  no  surprise  since  increases  in  P/d  "open 
up"  area  for  flow  through  the  crossflow  bank,  thereby 
decreasing  the  fraction  of  bypass  area.   The  increase 
in  R,  with  increasing  D.  is  less  than  that  of  P/d  , 
since  the  shell-to-outer-tube-limit  clearance  (Equation  3) 
also  increases  with  D.. 

Figure  12  is  a  graphical  representation  of  R, 

vs.  P/d  at  various  D- ,  in  the  same  manner  as  Figure  8. 
oi  ° 

It  should  be  noted  that  over  the  range  of  d  considered 

&      o 

(0.5  <  d„  <  1.5),  the  value  of  R,  will  decrease  about 
v       o      '  '  b 

5%  with  increasing  d  ,  a  much  larger  percentage  decrease 
than  can  be  predicted  for  J,  (and  subsequently,  heat 
transfer) . 
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3-  3.         3. 

Rb  "  W  l    ^P/do^  2  <V  3 


R 


e-c      Layout       a0       al       a2 


Square      .0277     -.0449    1.7242     .5932 
1  -  103     K2}*}**  .0542     -.0389    1.431      .4903 


Square 


Triangular    .0999     -.0324    1.1508     .392 


Square 

.0524 

-.0369 

1.4177 

.4878 

103  ■ 

■  105 

Rotated 
Square 

.0909 

-.032 

1.177 

.403 

Triangular 

.1504 

-.0267 

.9462 

.322 

TABLE  8   Correlation  Results  for  Bundle  Bypass  Factor,  R, 
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C.       CROSSFLOW  AREA  PRESSURE"  DROP 

Equation  A-23  gives  the  ideal  pressure  drop  for 
flow  across  the  tube  bank  between  baffles.   When 
multiplied  by  R,  and  R.,  the  effects  of  bundle  bypass 
and  baffle  leakage  are  accounted  for.   This  quantity 
has  been  correlated  in  the  same  manner  as  the  previous 
values.   The  friction  factor  correlations  of  Table  2 
were  used  along  with  other  corrections  for  tube 
orientation.   Table  9  is  the  result  of  these  overall 
correlations  for  the  two  flow  regimes  and  various 
orientations.   Note  that  the  exponent  xfi  of  the  basic 
equation  depends  on  flow  and  orientation  as  a  result 
of  the  friction  factor  correlation.   Hence,  the 
pressure  drop  can  not  be  expressed  only  as  a  function 
of  exchanger  parameters.   From  Table  9  it  is  obvious 
that  all  the  geometric  parameters  inversely  affect 
pressure  drop,  a  result  that  may  have  been  intuitively 
obvious,  but  the  extent  of  which  can  now  be  reasonably 
quantified.   Essentially,  the  increase  in  any  of  the 
five  basic  parameters  allows  for  more  area  for  fluid 
flow  with  less  restriction;  hence,  lower  pressure 
drops . 
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D.       WINDOW  AREA  PRESSURE  DROP 

Equation  A-23  gives  the  ideal  pressure  drop  for 
the  baffle  window  section  of  the  heat  exchanger.   When 
multiplied  by  R.  ,  the  effect  of  baffle  leakage  is  taken 
into  account.   Correlation  results  are  presented  in 
Table  10  for  differeing  tube  orientations.   Unlike 
crossflow  pressure  drop,  this  quantity  is  independent 
of  flow  regime. 

These  results  are  very  similar  to  the  crossflow 
pressure  drops  in  that  pressure  drop  is  inversely 
proportional  to  all  variables.   However,  in  the  window 
area  the  relative  importance  of  the  parameters  is  slightly 
different.   For  the  window  area,  shell  diameter  is  the 
predominant  influencing  factor  instead  of  pitch-to-diameter 
ratio.   The  relative  influence  of  baffle  cut  and  a  baffle 
spacing  ratios  is  nearly  the  same  for  both  area  pressure 
drops.   Tube  diameter  impact  is  relatively  the  smallest 
for  the  turbulent  crossflow  area  pressure  drop. 
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IV.   PERFORMANCE  AND  SIZING  OF  HEAT  EXCHANGERS 

A.       SENSITIVITY  ANALYSIS  OF  HEAT  TRANSFER  PER  UNIT 
PRESSURE  DROP 

For  any  given  application,  when  determining  the 
"optimal"  size  of  a  heat  exchanger,  some  feel  for  the 
overall  effect  of  parameter  variations  on  heat  exchanger 
performance  is  needed  in  order  to  judge  the  direction 
in  which  the  design  variables  should  be  altered  to  arrive 
at  a  reasonable  result.   The  performance  of  a  heat 
exchanger,  or  the  criterion  against  which  an  exchanger 
is  to  be  considered  optimal,  can  vary  widely  depending 
on  the  depth  of  the  analysis  that  is  being  performed. 
In  most  practical  applications  the  ultimate  criterion 
is  cost.   In  this  analysis,  cost  will  not  be  considered 
but  a  factor  which  can  influence  cost  to  some  extent 
will  be  analyzed. 

From  the  corrleation  results  for  heat  transfer 
coefficient,  Table  6,  and  the  basic  heat  balance: 


Q  -  AU0  AT,m  (8) 


consider  the  case  where  U   is  primarily  dependent  on  h  . 
Equation  (8)  then  can  be  written  as: 
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AT„   W  K- 
£m   s  1 


=  Ah 
o 


where  K,  is  defined  in  Table  6.   But  A  =  ud  NfcT/12  (ft2), 

1  0  tL 

and  the  results  from  Equation  (1)  (or  (2))  can  be  used 
for  N_  to  express  A  in  terms  of  d  ,  D- ,  P/d  .   Define: 

t  0    10 


rh    as 


iTlm  WsKl 


which  yields : 


7T  d  N   h 
o   t   o 


r./L  =  — u  L  u  (9) 

h  12 


This  quantity  is  a  non-dimensional  expression  for  heat 
transfer  per  unit  length  of  the  heat  exchanger.   From 
the  correlation  results  for  crossflow  and  window  area 
pressure  drops,  (Tables  9  and  10),  we  may  form, 
similarly: 

r  /L  =  --j-S (APb  ±    -    AP   .)/L  (10) 

s  /oc  ps 

This  is  a  simplified  expression  for  the  total  shell-side 
pressure  drop  per  unit  length  of  heat  exchanger.   Then 
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r, /r   can  be  expressed  as  a  quotient  of  those  correlations 

already  derived.   It  can  be  seen  from  all  the  correlation 

results  that  T-./T      is  not  only  dependent  on  the  geometry 

of  the  exchanger  but  flow  rate  also,  specifically, 

W  /u  .   For  this  analysis,  we  will  consider  only  turbulent 

flow  across  square-oriented  tube  bundles.   The  applicable 

correlation  results  were  selected,  the  ratio  T-./T 

'  h  p 

formed  and  that  ratio  was  graphed  as  a  function  of  the 

exchanger  parameters.   It  should  be  noted  that  a  Reynolds 

4 
number  of  10  was  used,  resulting  in  a  value  of 

tf  /y  =  2.026  x  10   at  the  mid-range  of  the  exchanger 

geometry  variables  (d  =  1.0,  P/d^  =  1.5,  £  /D.  =  .6, 

O  O  SI 

D°.  =  35).   Figures  13,  14,  15,  16  and  17  are  various  graphs 

of  I\/r  .   From  Figures  13-17  it  is  obvious  that  shell 
h  p  5 

diameter  D.  is  the  most  influential  parameter  affecting 
heat  exchanger  performance.   Figure  13  graphically 
illustrates  this  influence  for  several  values  of  baffle 
spacing  ratio,  I   /D. .   This  is  not  to  say  that  we  should 
always  optimize  an  exchanger  at  large  values  of  D-. 
Although  r,/r   increases,  this  graph  does  not  allow 
for  any  constraints  to  be  considered,  e.g.  maximum  shell- 
side  pressure  drop.   However,  we  may  conclude  that  where 
large  amounts  of  heat  are  to  be  transferred  with  low 
shell-side  pressure  drops  as  a  requirement,  "large"  shell 
diameters  will  be  required. 
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Figure  14  illustrates  the  influence  of  baffle 

spacing  on  T./T      at  several  values  of  D. .   This  parameter 

r  °  h  p  i 

also  favorably  affects  the  heat  exchanger  performance 
but  certainly  not  to  the  extent  of  shell  diameter.   Again, 
the  realization  of  the  absence  of  a  pressure  drop  con- 
straint must  be  accounted  for. 

Figure  15  illustrates  the  relationship  between 
I\  /r   and  tube  diameter.   This  is  the  only  parameter  (of 
those  considered)  which  adversely  affects  I\  /Y    .   The 
most  basic  explanation  for  this  effect  is  that  increases 
in  tube  diameter  will  more  quickly  decrease  heat  transfer 
than  they  will  pressure  drop. 

Figure  16  and  17  relate  the  effects  of  baffle 

cut  ratio,  I   /D-  and  pitch-to-diameter  ratio,  P/d  ,  to 
c   1  o 

IV /r  .   It  is  obvious  that  at  the  conditions  analyzed 
h  p  J 

(turbulent  flow,  square  tube  orientation) ,  neither  para- 
meter has  any  significant  effect  on  heat  exchanger 
performance. 

The  five  figures  shown  can  be  of  assistance 
when  determining  exchanger  size  for  a  given  "job".   The 
qualitative  effect  (or  lack  thereof)  of  the  five  para- 
meters will  help  guard  against  "false  starts"  in  any 
iterative  sizing  scheme. 
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FIGURE  16    r,  /r   vs  2  /D. 
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log-rt  r,  /T  Turbulent  Flow  Square  Orientation 

s10   h   p 
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0.5  W                                   - 
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FIGURE    17      rh/r      vs    P/dQ 
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B.       PROPOSED  SIZING  ROUTINE  FOR  CASE  OF  SHELL -SIDE 
PROPERTIES  ONLY 

For  the  case  where  the  shell-side  heat  transfer 

coefficient  is  the  controlling  factor  (the  tube  side 

coefficient  has  little  impact  on  overall  U  ) ,  it  is 

desirable  to  determine  the  exchanger  geometry  for  a 

given  "job":   total  Q  transferred; AP   specified;  W  , 

AT   (shell-side),  and  fluid  properties  all  given.   We 

may  desire  to  "optimize"  the  heat  exchanger  on  minimum 

volume  where  volume  is  simply  expressed  as: 

it  D.2 

V  =  i-  L 


From  the  basic  heat  balance  equation, 


tt  d  N.  L 
Q  =  AUo  AT.  .   U   =  h  .    A  =  2 — I —      (ii) 

x       0     £Hl     O      O  ,  ~ 


yields 


7r  d  N  L 

Q  =  AT£m     12    ho  ^12^ 


With  knowledge  of  the  fluid  properties  we  may  form  the 
quantity  I\  : 

Q 

rh  =  (13) 

n    AT,   K^ 
Jim   1  s 


65 


where  Ki  is  specified  in  Table  6,  correlation  results 
for  shell-side  coefficient.   Then 


.  7T  d  N.  L 

rh/L  =  °— ^—  h  (14) 

n         12     ° 


We  must  now  begin  to  specify  the  geometry  of  the 
exchanger.   First,  pick  d  ,  P/d  ,  D.  as  the  beginning 
points.   From  equation  A- 12  for  Reynolds  number:- 


•g 


12d  W 
R   =  2_§. 

s  m 


(15) 


where 


W         a        a 

—  c   (d  )  1    (P/d  )  2  (£C/D.)~3  (D-)"4 
0    0     K         0J  s   ly     **  i' 

ys 


c   =  62.224,   a,  =  1.037,   a-  =  -  1.915, 
o  1  z 


a.  =  1.0,   a.  =  -  1.933 
3  4 


we  may  either  determine  &  /D-  (if  we  desire  to  operate 

at  a  certain  R  )  or  pick  £  /D.  and  determine  the  flow 
e„^    r     si 
s 

regime  (laminar/transition  or  turbulent) . 

Once  this  condition  has  been  finalized,  the 
proper  correlation  for  h   is  selected  from  Table  6, 
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knowing  the  selected  tube  orientation.   The  last  para- 
meter to  be  determined  is  I   /D . .   As  Bell  has  stated, 
for  "well-designed"  exchangers,  J   should  be  about  1.0, 
which,  from  Figure  6,  yields  a  value  of  I   /D.    of  about 

*—     J- 

0.25.   The  initial  value  of  h  may  now  be  calculated 

from  the  correlation  results  of  Table  6  and  exchanger 

length  L  determined  from  Equation  10,  using  the  proper 

N  correlation  for  tube  orientation.   Volume  may  then 

be  easily  calculated. 

Pressure  drop  calculations  can  be  made  from  the 

correlation  results  of  Tables  10  and  11,  knowing  the 

exchanger  length  L.   This  must  then  be  checked  against 

the  allowed  AP  .   If  AP     .  <  AP        ,  we  may 

calc      allowed 
decrease  the  size  of  the  exchanger  to  utilize  more  of 

the  allowed  AP  .   This  is  most  easily  accomplished  by 

decreasing  D-  and  recalculating  h  ,  L,  AP   and  V,  to 

arrive  at  some  minimum  volume. 

If  AP,.   ,   >  APc-  ^T    ,,  the  procedure  is 
scalc     Allowed*     r 

reversed,  i.e.,  increase  D-  and  recalculate  the 
necessary  quantities  until  a  convergence  within  some 
desired  bound  is  obtained.   As  in  the  case  in  any 
sizing  procedure,  some  notes  of  caution  and  guidelines 
are  in  order.   When  increasing  shell  diameter  (or  any 
of  the  other  parameters) ,  caution  must  be  taken  that  the 
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flow  regime  remains  as  desired.   It  is  also  possible  to 
"fine  tune"  the  exchanger  by  altering  the  other  geometric 
values.   If  the  difference  APS(;alc  -  APsallowed  is  small, 
we  may  desire  to  change  the  baffle  cut  or  baffle  spacing 
ratios.   Small  changes  in  these  parameters  will  also 
result  in  changes  to  the  shell  side  pressure  drop  and 
heat  transfer  coefficient,  the  effect  of  which  can  be 
examined  by  the  exponents  of  the  applicable  correlation 
result.   We  may  also  use  Bell's  insight  into  the  heat 
transfer  corrections  to  assure  that  the  heat  exchanger 
is  reasonably  well  designed.   For  the  heat  transfer 
correction  factors  J,  ,  J  ,  J.   Bell  gives  the  following 

ranges  as  good  guidelines: 


J   *   1.0 
c 


0.7   <   J,   <   0.9 

D 


0.7   <   J    <   0.8 


From  Figure  9,  10,  and  11,  we  can  pick  reasonable  values 

of  S.   /D.  and  £  /D.  as  starting  points  for  the  exchanger 

c   l      s   i  °  r  ° 

geometry  and,  should  changes  be  made,  ensure  that  we 
remain  in  a  normally  expected  range  for  those  values. 
The  use  of  those  figures  as  aids  to  sizing  should  prevent 
haphazard  variations  and  point  to  a  reasonable  solution 
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FIGURE  18   Simplified  Logic  Sequence  for 
Proposed  Sizing  Routine 
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in  an  expedient  manner.   Figure  13  gives  a  simplified 
flow  chart  of  the  logical  sequence  for  this  proposed 
sizing  routine. 

C.       EXAMPLE  OF  PROPOSED  SIZING  ROUTINE'  WITH 
COMPARISON  TO  THE  DELAWARE  METHOD 

The  preceding  sections  of  this  analysis  have  dealt 
with  the  actual  results  of  the  non-linear  regression  used 
to  simplify  the  heat  transfer  coefficient  and  pressure 
drop  calculations.   In  this  section,  a  hypothetical 
example  for  determining  exchanger  size  is  given  to  show 
the  actual  mechanics  of  the  proposed  sizing  routine  and 
to  compare  those  results  with  the  calculations  of  the 
Delaware  Method.   This  comparison  is  done  for  the  case 
where  only  the  shell-side  heat  transfer  coefficient  and 
pressure  drop  are  of  concern  to  the  designer. 

The  "job"  considered  in  this  example  is  as 
follows : 

1)  shell-side  mass  flow  rate  of  water: 

W  =  250,000  lbm/hr 
s      ' 

2)  heat  water  from  50°F  to  90°F 

3)  maximum  shell-side  pressure  drop: 
APS  <_  8.0  psi 
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The  following  fluid  properties  are  assumed: 


1) 
2) 

3) 
4) 
5) 


1.0  btu/lbm  -  °F 


p    =1.65  lbm/hr-ft  (at  an  assumed 


T   n  =  100°F) 
wall        ' 


ys  =  2.445  lbm/hr-ft  (at  Taye  =  70UF) 


Prs  =  7.05 

ps  =  62.075  lbm/ft 


We  desire  to  operate  in  the  turbulent  flow  regime 

and  must  check  the  Reynolds  number  for  our  exchanger. 

Picking  tube  diameter  d  =0.75  in,  pitch-to-diameter 

ratio  P/d  =  1.33,  baffle  spacing  ratio  £  /D.  =  0.75 
o  r     ©         si 

and  shell  diameter  D.  =  25.0  in,  a  check  of  Reynolds 
number  using  Equation  (15)  yields  Re   =  7229,  which 
ensures  turbulence.   Using  a  square  tube  orientation, 
baffle  cut  ratio  I   /D.  =  0.25,  we  select  the  applicable 
correlation  result  from  Table  6: 


-°  =  W 


fyf   "\  "•  353 


t  s 


.212 


(4.07)(dQ) 


.015 


(P/dQ) 


.4492 


.376 


U  /DO 

^  c   1 


.  cvv 


1.032 


(j>0 


where 


K-,  =  c   (P   ) 

1    s  *■  r  J 


2/3 


^0.14 


u 
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Using  the  values  given  results  in  a  heat  transfer 
coefficient  of: 


h   =  463.5  Btu/hr-ft3-°F 
o 


From  Q  =  AU   AT „   =  Ah  AT„  ,  where 
x     o    in  o   )ln' 


2.17526 


.  7T  d  Nr  L  D. 

A  = — - —  and  N   =  0.3559  — » 

12  Z  d  " 


The  length  L  of  the  heat  exchanger  is  calculated 


L  =  6.99  ft, 
then  volume 

V  =  23.83  ft3  . 

Using  the  correlation  results  for  crossflow  and  window 
area  pressure  drops  from  Tables  9  and  10  respectively, 
yields  the  total  shell-side  pressure  drop: 


AP  =  .29  psi 
s 


As  a  comparison,  utilizing  the  Delaware  Method  (Reference 
1)  ,  the  following  is  presented: 
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Proposed  Method  Delaware  Method 

h   (Btu/hr-ft2-°F)         463.5  454.8 

o  K        ' 

L  (ft)  6.99  7.05 

V  (ft3)  23.83  24.03 

AP  . (psi)  0.29  0.25 


Although  the  results  are  quite  close,  we  are  far  from 
"optimized"  in  using  the  allowed  shell-side  pressure  drop 
We  may  decrease  shell  diameter  D.  to  most  easily 
accomplished  this.   For  the  next  interation,  pick 

D.  =  15.25  in.   This  "odd"  size  is  picked  only  because 

i 

Bell's  tables  in  Reference  (1)  list  all  the  applicable 
parameters  for  this  size  (D  t«,  N  ,  etc.)  and  inter- 
polation between  those  listed  values  would  not  result  in 
a  meaningful  comparison. 

A  recalculation  of  the  heat  exchanger  using  this 
new  value  of  D.  yields  the  following  result: 

Proposed  Method      Delaware  Method 

h  785.16  789.0 

o 

L  12.09  11.69 

V  15.34  14.83 

AP_  3.62  3.61 


73' 


These  results  are  as  agreeable  as  those  of  the 
first  iteration,  yet  the  exchanger  may  be  further  de- 
creased in  size.   A  third  iteration,  using  D-  =  13.25  in 
yields : 

Proposed  Method      Delaware  Method 

h  907.8  917.5 

o 

L  14.20  14.61 

V  13.6  13.99 

AP  7.53  7.27 

s 

These  results  show  that  we  have  come  close  to  a 
preliminary  size  for  the  heat  exchanger  that  meets  the 
required  "job".   Some  comments  concerning  the  previous 
calculations  are  noteworthy. 

In  decreasing  shell  diameter  in  these  iterations, 
the  recalculation  of  h  by  the  proposed  method  is 
accomplished  simply  by  multiplying  the  previous  value 
by  the  ratio  of  shell  diameters  raised  to  the  applicable 
correlation  exponent.   Pressure  drop  calculations  are 
estimated  in  much  the  same  manner,  except  that  the 
crossflow  and  window  area  terms  must  be  done  separately. 
The  only  "lengthy"  calculation  involves  that  of  exchanger 
length,  L. 
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For  the  sizing  determination  by  the  Delaware 
Method,  no  less  than  twenty-one  (21)  separate  calculations 
and  seven  (7)  graphical  readings  must  be  done  for  each 
iteration.   For  this  proposed  method  only  five  (5) 
calculations  were  necessary:   heat  transfer  coefficient 
(Table  6) ,  overall  length  (Equation  12) ,  pressure  drops 
for  crossflow  and  window  sections  (Tables  9  and  10) , 
and  total  pressure  drop  (Equation  A-25). 

This  quick,  handy  method  is,  in  the  opinion  of 
this  author,  much  more  preferable  to  the  original  method 
for  preliminary  heat  exchanger  sizing  and  "back-of-the- 
envelope"  calculations.   It  is  realized  that  no  simplified 
heat  exahanger  model  can  replace  the  more  sophisticated 
computer  optimization  methods  that  have  been  and  are 
being  developed.   However,  in  a  relatively  short  time, 
the  designer  can  achieve  a  realistic  estimate  of  the 
heat  exchanger  performance  and  size,  and  can  use  these 
results  as  possible  inputs  to  the  more  elaborate  routines. 
The  closer  an  initial  estimate  is  to  the  actual  optimum, 
the  more  time  and  money  can  be  saved  in  the  more 
sophisticated  (and  expensive)  methods.   Minimization 
of  cost  will  almost  always  be  a  major  goal. 
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APPENDIX  A 

Basic  equations  used  for  heat  transfer  coefficient 
and  pressure  drop  calculations. 
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(8) 


S   ,     =    tt/2    D.    6   ,     [1    -    3/2] 
sb  l     sb 
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(16)  J      =    0.55    +    0.72    F 

*■  c  c 


-2.2r 


(17)  Jz   =   0.44    (1    -    rs)    +    [1    -    0.44(1    -    rs)]e 


im 


(18)  J,     =   e 
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v      J  o  ideal      c      £      d 
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(24)  AP 


w,i 
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s       2 cw 

2    2      p      S      S 
5c      s      m     w 


(25)  APS    =    C(Nb    +    1)    APb>i    Rb    ♦    Nb    4P„   j^ 
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